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We present a non-isothermal mesoscopic model for investi¬ 
gation of the phase transition dynamics of thermorespon¬ 
sive polymers. Since this model conserves energy in the 
simulations, it is able to correctly capture not only the 
transient behavior of polymer precipitation from solvent, 
but also the energy variation associated with the phase 
transition process. Simulations provide dynamic details of 
the thermally induced phase transition and confirm two 
different mechanisms dominating the phase transition dy¬ 
namics. A shift of endothermic peak with concentration is 
observed and the underlying mechanism is explored. 

Thermoresponsive polymers (TRPs) have attracted increasing 
attention in the last two decades because of their great po¬ 
tential applications in various chemical and biological sys- 
temslffil, i.e. , controlled drug delivery, smart materials, biosep¬ 
arations and filtration. Most applications of TRP have re¬ 
lied on a drastic and discontinuous c hang e of their solubil¬ 
ity in given solvents with temperatureCEl. In particular, the 
temperature-composition diagram of TRP involves a misci¬ 
bility gap. Depending on the miscibility gap if it appears at 
low or high temperatures, the critical temperature Tc is known 
as the lower critical solution temperature (LCST) or the up¬ 
per critical solution temperature (UCST), respectively. LCST- 
type TRPs are hydrophilic and highly mixed with the sur¬ 
rounding solvent at low temperatures, but become hydropho¬ 
bic and precipitate from the solvent above the critical phase 
transition temperature, while UCST-type TRPs exhibit the op¬ 
posite behavioi®. The underlying mechanism of this solubil¬ 
ity transition with temperature is related to the role of hydro¬ 
gen bonds®. For LCST-type TRP at low temperature T <Tc, 
hydrogen bonds are generated between solvent and polymer 
molecules. Therefore, the polymers show hydrophilic prop¬ 
erties and can be easily dissolved into the solvent. However, 
when the temperature is increased above the critical tempera¬ 
ture T > Tc, those solvent-polymer hydrogen bonds are dis¬ 
turbed and polymer-polymer hydrogen bonds dominate the 
dynamics, which makes the polymer become hydrophobic and 
precipitate from the solvent. 

In practical applications of TRP, temperature-sensitive mi¬ 
crogels/micelles are often used for the functional element®. 
The major building block for these temperature-sensitive mi¬ 
crogels is TRP. Among them, poly(N-Isopropylacrylamide) 
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(PNIPAM) is the most investigated material and was exten¬ 
sively used for the construction of temperature-sensitive mi¬ 
crogels. Specifically, PNIPAM has a LCST around 32°C be¬ 
tween room and body temperatures, which makes it a promi¬ 
nent candidate in biomedical applications®. As a matter of 
fact, the applications of TRP highly depend on the evolution of 
the microstructure of microgels in the phase transition process. 
Predicting the performance of TRP-based materials requires 
a deep understanding of the thermally induced phase transi¬ 
tion dynamics. Usually, experiments are able to observe the 
coil-to-globule transition of polymers by measuring the light 
transmittance rate, and study the static microstructure of mi¬ 
crogels using NMR, light scattering and transmission electron 
microscopy®. It is well-known that the static microstructure of 
LCST-type microgels is swollen at low temperatures and col¬ 
lapsed at high temperatures, as shown in Fig. To this end, 
some theories have been developed for understanding the ex¬ 
perimental observations at the molecular level®. However, it 
is difficult for both experiment and theory to provide dynamic 
details of the transition process, which is very important for 
clarifying the phase transition of TRP. Alternatively, compu¬ 
tational simulation techniques are able to provide details of 
these dynamic processes and can be used to assist in silico 
design for specific applications of TRP. 

The typical diameter of a single microgel bead ranges be¬ 
tween 50 nm and 5 jurr^, and the time scale for the phase 
transition of individual PNIPAM-based microgel beads is in 
the order of 100 n^. Simulating a dynamic process lasting 
for hundreds of nanoseconds with a length scale of several 
microns is difficult for conventional atomistic methods, but 
it is in the comfortable temporal and spacial scales of meso¬ 
scopic approachesESI. Therefore, in the present work, we de¬ 
velop a non-isothermal mesoscopic model based on dissipa- 



Fig. 1 A single temperature-sensitive microgel bead consisting of 
LCST-type TRP chains showing (a) a fully swollen state at low tem¬ 
perature T <Tc, and (b) a fully collapsed state at high temperature 
T>Tc. Various polymer chains are visualized with different colors. 





live particle dynamics (DPD) to investigate the thermally in¬ 
duced phase transition of TRR DPD is a particle-based meso¬ 
scopic approach, which is usually considered as a coarse¬ 
grained molecular dynamics (MD) model^. The classic DPD 
method was designed for simulating isothermal hydrodynam¬ 
ics, which is not valid for non-isothermal processes because 
of the violation of energy conservation^l. To conserve the en¬ 
ergy of the system, an extension of DPD was developed by 
including the mesoscopic energy equationEESl. The energy- 
conserving DPD model is known in the literature as eDPD, 
and it has been demonstrated that eDPD conserves the energy 
of fluid systems in simulations and can capture the correct 
temperature-dependent properties of fluids^El. In this paper, 
we extend the eDPD framework to modeling the temperature 
sensitivity of TRP (for details on the eDPD formulations, see 
ESI^). 

The potential between TRP and solvent is sensitive to the 
temperature changes, and the polymer-solvent interaction pa¬ 
rameter X is a function of temperaturel^. In the DPD method, 
the Flory-Huggins y-parameter is linear with respect to the 
excess repulsion which is defined by Aa = asp — ass 

where 5' represents solvent and p stands for polymer. To 
model the thermally induced phase transition of TRP, we de¬ 
fine the excess repulsion Aa as a function of temperature to 
consider the temperature-dependence of the Flory-Huggins y- 
parameter. In practice, we take the repulsion parameters be¬ 
tween particles as ass{T) = app{T) = ISksT /p, and the cross 
terms asp{T) =15kBT /p+Ao + AA/[1.0-l-^vp(— x- {T — Tc))] 
containing a sharp change by AA at T = (see Fig. S2, 
ESIt). Since the conservative force between particles is given 
by Ffj = aij{T){l — rij/rc)eij and the corresponding poten¬ 
tial is Uij = ^aij{T)rc{l — rijjrc)^, the pair potential between 
particles changes with temperature because of the variation of 
repulsive coefficient atj. To satisfy the conservation of energy, 
the change of potential energy is considered to be balanced by 
a change of internal energy. Specifically, the total energy for 
each pair is considered invariable and its variation is zero upon 
time integration, i.e., AEij = AUij -f Aci — Qi + Aej — Qj = 0, 
where Aet = CyATi and Aej = CyATj are the changes of inter¬ 
nal energy of particles i and 7 , and Qt and Qj represent their 
net heat fluxes. 

We consider an eDPD system containing one microgel bead 
in solution, and take LCST-type TRP as an example. Here, the 
results are interpreted in terms of the reduced DPD units, un¬ 
less specified otherwise. For applications to specific materials, 
we refer interested readers to Refs. [T3l[T4]for parameteriza¬ 
tion of DPD systems. The microgel bead is made up of many 
cross-linked linear polymer chains (see Fig. SI, ESIf). Each 
polymer chain consists of 50 eDPD particles sequentially con¬ 
nected by harmonic springs, and cross-links with a density of 
approximately 3% total bonds are randomly distributed in the 
microgel bead. The eDPD system of a TRP microgel bead 
surrounded by solvent particles is initialized at a low temper¬ 


ature To = O.STc. Simulations involving half a million parti¬ 
cles are performed using a GPU-accelerated DPD user^'^^0 
package!^. Since the particle system is constructed with ran¬ 
dom initial configurations, we run the eDPD simulations at 
To = O.STc for 100 time units to obtain the thermal equilibrium 
state. Then, similarly to the differential scanning calorimetry 
(DSC) experiments, the temperature of the eDPD system is 
increased linearly as a function of time, i.e., from Tq = O.ST^ 
to 7i = lATc within 1500 time units. 

An obvious observation in the eDPD simulation is the con¬ 
figurational change and associated change in size of the TRP 
microgel bead. To quantify the deformation of the microgel 
bead during the heating process, we computed its instanta¬ 
neous gyration radius Rg. Fig. [^a) shows the time evolu¬ 
tion of Rg during heating from O.ST^ to lATc, which con¬ 
tains a significant decrease of Rg corresponding to a phase 
transition between T = Tc and T = lATc. More specifically, 
eight configurations of a large microgel bead (5 wt%) along 
the Rg curve are presented in Figs, [^bi-bg), and four snap¬ 
shots for a small microgel bead (1 wt%) are shown in Figs. 
[^Ci-C 4 ). The corresponding points are marked with symbols 
on the curve of Rg shown in Fig. |^a). At temperatures be¬ 
low the critical temperature of phase transition, i.e., T < 
the TRPs are hydrophilic leading to a fully swollen state of 
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Fig. 2 (a) Evolution of gyration radius Rg of LCST-type thermore- 
sponsive microgels with concentrations of 1 wt% and 5 wt% during 
heating, and (C 1 -C 4 ) show cross sections of their transient 

microstructure corresponding to the changes of Rg. 


If] 


2 | 


2 





































the microgels, as displayed in Figs. Hbi, Cl). A swollen 
microgel bead has the maximum volume and corresponds to 
the maximum gyration radius. Thus, the curves of Rg have a 
plateau when T <Tc. However, as the temperature increases, 
the hydrophilic-hydrophobic transition occurs near the critical 
temperature T ^Tc, above which the TRPs become hydropho¬ 
bic, and hence the microgel bead starts to collapse until it turns 
into a compact collapsed globule as shown in Figs, [^bg, C 4 ). 

For self-aggregation of TRP in the coil-to-globule phase 
transition process, two different mechanisms dependent on the 
size of TRP molecules dominate the dynamicsSS. Small TRP 
molecules undergo an “all-or-none” process while large TRP 
molecules behave as if they consist of quasi-independent “do¬ 
mains”® Figs. ibi -bg) and (C 1 -C 4 ) provide dynamic details 
of the phase transition process. Our simulations confirm that 
a small microgel bead (1 wt%) has the “all-or-none” coil-to- 
globule transition, which is a relatively simple process (Fig. 
[^ci-C 4 )). However, a larger microgel bead (5 wt%) has many 
“independent domains” that start their self-aggregation pro¬ 
cesses simultaneously. In particular, porous structures can be 
observed at the beginning of coil-to-globule transition of the 
microgel bead, as shown in Fig. [^b 2 ). Those pores may trap 
solvents and then merge them inside the microgel bead (Fig. 
I3b4 -b6)). In general, a big solvent droplet trapped inside a hy¬ 
drophobic microgel bead is unstable. The droplet will finally 
escape, and then a compact collapsed globule is observed in 
Fig.|3bg). 

Experiments^^® reported that LCST-type TRP undergo a 
coil-to-globule transition upon temperature increase, which 
is accompanied by cooperative heat absorption and results in 
an endothermic peak in the differential scanning calorimetry 
(DSC) thermogram. In particular, the endothermic peak de¬ 
pends on polymer concentration and shifts slightly to lower 
temperatures with concentrating polymer solutions®. Figure 

shows the thermal events that require energy flow to eDPD 
systems with a fixed scanning rate (i.e., 4.0 x lO^'^Tc per time 
unit) for different concentrations (1 wt%, 2 wt%, 5 wt%, 
10 wt%). We observe a shift of the endothermic peak to lower 
temperatures as the concentration increases. Since the effects 



Fig. 3 Thermogram analog to differential scanning calorimetry 
(DSC) heating curve with a fixed scanning rate at different concen¬ 
trations (1 wt%, 2 wt%, 5 wt%, 10 wt%). 


of ionic additives and PH are excluded in the present model, 
the only reason responsible for the shift is the difference of 
self-aggregation dynamics. At the beginning of the phase tran¬ 
sition, the nucleation dynamics of TRP is reversible. For small 
microgels, the “all-or-none” process needs more time to ini¬ 
tialize the phase transition. However, many “independent do¬ 
mains” of larger microgels start their nucleation processes si¬ 
multaneously, which reduces the delay of phase transition. As 
a demonstration case of the proposed model, our simulations 
confirm that small TRP microgels in solution experience an 
“all-or-none” process while large TRP microgels are domi¬ 
nated by a “domain” mechanism®(the independent formation 
of cooperative units). The different dynamic processes con¬ 
tribute to the shift of endothermic peak in DSC thermograms. 

In summary, we present a non-isothermal mesoscopic 
model that can be used to simulate thermally induced phase 
transition of TRP and to provide dynamic details of the phase 
transition process. Since the model conserves the energy of 
system in the simulations, it is also able to capture the un¬ 
derlying energy variation associated with the phase transition. 
This mesoscopic model is a promising candidate for modeling 
thermally induced phase transition of various thermorespon- 
sive polymers and can assist in silico design for engineering 
and biomedical applications of thermoresponsive materials. 
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